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"Method and apparatus for magnetic field detection" 
Technical Field 

The present invention relates to a method and apparatus for 
measurement of a magnetic field, and in particular relates to rotating a 
5 magnetic sensor relative to the magnetic field to be measured. 

Background Art 

A gradiometer is a device for measuring the difference between two 
signals, with an aim to facilitate rejection of common mode noise signals and 
io improved reduction in errors that arise from misorientation of the measuring 
instrument. 

The ability of a gradiometer to measure small magnetic field gradients in 
the presence of a large background magnetic field is determined by its ability to 
reject large common mode signals, and the intrinsic noise and drift 

15 performance of the field detectors employed, The first of these characteristics 
is called the common mode rejection ratio (CMRR). In conventional 
gradiometers the CMRR is determined by the precision with which the field 
sensors comprising the gradiometer are balanced. That is, the CMRR is 
determined by the relative sensitivities of the field sensors, and the accuracy to 

20 which the field sensors can be positioned in precisely parallel planes. This 
fundamental nexus between the precision with which conventional 
gradiometers can be manufactured and the ultimate CMRR performance has 
until now limited the CMRR to figures of the order of 10" 4 in practical 
instruments. 

25 Further, the DC and low frequency accuracy of state of the art 

gradiometers is limited by the DC and drift performance of the field detectors 
employed, and is also limited by the low frequency noise performance of the 
field detectors. Furthermore, the achievable noise performance can be 
dependent upon the magnitude of the background homogeneous field; being 

30 determined by microphonics which arise from vibrations causing randomly 
varying misalignment of the axes of symmetry of the two field sensors. For 
example, in the case of magnetic gradiometry two appropriate magnetic 
sensors are the flux gate and the superconducting quantum interference device 
(SQUID). Flux gates have considerable drift and noise in comparison to 

35 SQUIDs, but state of the art SQUID systems are incapable of providing data 
about the absolute value of the magnetic field or gradient. 
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The inability of state of the art SQUIDs to measure the absolute value of 
the magnetic field or gradient arises for the following reasons. The radio 
frequency (rf) SQUID is essentially a ring of superconducting material, the ring 
being interrupted by a Josephson Junction. When the superconducting ring is 

5 energised by an inductively coupled resonant rf-oscillator, tunnelling of 
electrons takes place at the junction and a periodic signal, being a function of 
flux through the ring, can be detected across the junction. The periodic signal 
is substantially a triangular waveform, usually having a period (AB) in the order 
of a nanotesla. Therefore, in order to yield a sensitivity in the femtotesla range, 

10 the SQUID is operated in a nulling bridge mode, or flux locked loop (FLL) 
mode. In this mode, magnetic flux is fed back to the SQUID so as to cause the 
output voltage to remain relatively constant. The feedback voltage, being 
proportional to the difference between the applied flux and the quiescent flux 
level, gives a highly accurate measurement of relative magnetic flux. The 

15 feedback voltage V can therefore be written as 

V = M(A eff B + u) (1) 

where 

M is a constant in a specific SQUID system; 
A e ff is the effective area of the SQUID; 
20 S is the applied magnetic field; and 

u is the quiescent flux. 

However, the quiescent flux u is unknown. Thus, conventional rf- 
SQUIDs provide only relative measurements of magnetic field, up to an offset 
(which can usually be held constant for the duration of the measurement 

25 procedure) which is an integer multiple of a field value determined by the 
sensitivity of the SQUID and the fundamental flux quantum. The SQUID, 
therefore, is not an absolute field detector and state of the art gradiometers 
cannot provide absolute value measurements of the gradient components. 

A further problem arises in tensor gradiometers. The gradient of a field 

30 in 3 dimensions is a second rank tensor, and accordingly has 9 components. 
In the case of magnetic fields Maxwell's equations impose restrictions which 
imply that only 5 of the gradient components are independent. To determine 
the gradient tensor of a magnetic field it is therefore necessary to measure all 5 
of these components. State of the art tensor gradiometers achieve this by 

35 deploying a suite of at least 5 field sensors or planar gradiometers, or a 
combination of field sensors and planar gradiometers. 
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However, a problem which arises from the use of a plurality of field 
sensors or gradiometers to measure the gradient components, is that the 
sensors or gradiometers distort the field in their neighbourhood. Consequently, 
each sensor or gradiometer adversely affects the measurement of the field by 
5 the other sensors or gradiometers which are deployed to measure other 
gradient components. Such a "cross-talk" effect limits the accuracy of the 
gradient measurement. 

Any discussion of documents, acts, materials, devices, articles or the like 
which has been included in the present specification is solely for the purpose of 
io providing a context for the present invention. It is not to be taken as an 
admission that any or all of these matters form part of the prior art base or were 
common general knowledge in the field relevant to the present invention as it 
existed before the priority date of each claim of this application. 

Throughout this specification the word "comprise", or variations such as 
15 "comprises" or "comprising", will be understood to imply the inclusion of a 
stated element, integer or step, or group of elements, integers or steps, but not 
the exclusion of any other element, integer or step, or group of elements, 
integers or steps. 

Throughout this specification, the terms 'superconducting material', 
20 'superconducting device' and the like are used to refer to a material or device 
which, in a certain state and at a certain temperature, is capable of exhibiting 
superconductivity. The use of such terms does not imply that the material or 
device exhibits superconductivity in all states or at all temperatures. 

25 Summary of the Invention 

According to a first aspect, the present invention provides a method of 
field detection comprising: 

providing a gradiometer in a field, the gradiometer having at least first 
and second field vector sensors connected in a differencing arrangement; and 
30 controllably altering a position of the at least first and second field vector 

sensors relative to the field during operation of the gradiometer. 

According to a second aspect the present invention provides a field 
detection device comprising: 

a gradiometer, the gradiometer having at least first and second field 
35 vector sensors connected in a differencing arrangement; and 
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means for controllably altering the position of the at least first and second 
field vector sensors relative to a field during operation of the gradiometer. 

In preferred embodiments of the present invention, the at least first and 
second field vector sensors of the gradiometer are substantially axially aligned, 

5 thus forming an axial gradiometer. That is, the at least first and second field 
vector sensors are preferably positioned relative to each other such that first 
and second field vectors sensed by the respective first and second field vector 
sensors are substantially co-axial. However, the present invention may also be 
applied in respect of other devices such as planar gradiometers. 

io In particularly preferred embodiments of the present invention, the 

position of the at least first and second field vector sensors is controllably 
altered by rotating the at least first and second field vectors about an axis of 
rotation. Such rotation may be continuous, as in the case of a motor driving the 
gradiometer at a substantially constant revolutional velocity about the axis of 

15 rotation. Alternatively, the rotation may be piecewise about the axis of rotation, 
as in the case of manual rotation of the gradiometer through a given angle 
about the axis of rotation. Measurements may be obtained from the 
gradiometer before, during, and/or after such a manual rotation. 

In preferred embodiments of the invention, the axis of rotation is 

20 positioned substantially perpendicular to the co-axial first and second field 
vectors. The axis of rotation preferably passes substantially between the first 
and second field vector sensors, and substantially equidistant from the first and 
second field vector sensors. 

In accordance with preferred embodiments of the present invention, the 

25 gradiometer is made to rotate so that the direction of the sensitive direction of 
the field vector sensors or their sensitivities is changed in a known way. 

In accordance with preferred embodiments of the invention, the 
gradiometer may be an axial magnetic gradiometer, and may be formed from a 
pair of SQUIDs, flux gates or superconducting pick-up loops with their 

30 sensitivity vectors lying substantially in a nominal x«y plane. In preferred 
embodiments, such an axial gradiometer is rotated about a nominal z-axis to 
form the rotating gradiometer. 

While the use of an axial gradiometer may be favoured in 
implementations of the present invention because of its intrinsic noise 

35 immunity, the advantages of rotating a gradiometer are not limited to the axial 
configuration. 
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It has been realised by the present inventors that a rotating axial 
magnetic gradiometer returns information about the value of the gxy component 
and linear combinations of the values of the g,* and g yy components of a field 
gradient tensor, in addition to information about the B x and B y field components. 
5 In preferred embodiments of the invention, such information is retrieved by 
transforming the voltage waveform from the field vector sensors into the Fourier 
domain. 

Furthermore, it has been realised that it is possible to obtain the 
complete magnetic gradient tensor by deploying three such rotating 
10 gfadiometers, configured in such a way that their respective axes of rotation 
are not parallel. 

Accordingly, in a third aspect, the present invention provides a method of 
obtaining a complete magnetic gradient tensor of a magnetic field, the method 
comprising: 

15 providing at least three axial gradiometers in the magnetic field such that 

an axis of each axial gradiometer is not parallel to an axis of any other one of 
the at least three gradiometers; and 

controllably altering a position of each of the at least three gradiometers 
relative to the magnetic field during operation of the gradiometer. 

20 Preferably the step of controllably altering comprises rotating each of the 

at least three axial gradiometers about a respective axis of rotation of that 
gradiometer. The axis of rotation of each gradiometer is preferably not parallel 
to an axis of rotation of any other one of the at least three gradiometers. 

Similarly, in a fourth aspect, the present invention provides a device for 

25 obtaining a complete magnetic gradient tensor of a magnetic field, the device 
comprising: 

at least three axial gradiometers positioned such that an axis of each 
axial gradiometer is not parallel to an axis of any other one of the at least three 
gradiometers; 

30 means for controllably altering a position of each of the at least three 

gradiometers relative to the magnetic field during operation of the gradiometer. 

Preferably the means for controllably altering comprises means for 
rotating each of the at least three axial gradiometers about a respective axis of 
rotation of that gradiometer. The axis of rotation of each gradiometer is 

35 preferably not parallel to an axis of rotation of any other one of the at least 
three gradiometers. 
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In preferred embodiments of the third and fourth aspects of the invention, 
the at least three gradiometers are rotated at a substantially constant velocity 
about their respective axes of rotation. Such embodiments produce data 
suitable for conversion to the Fourier domain. 

5 Embodiments of the third and fourth aspects of the invention effectively 

provide a single measurement of the true value of each of the three off- 
diagonal components of the gradient tensor, three measurements of the 
diagonal components and two measurements of each of the field components. 
Because more than the minimum of five tensor components are obtained, the 

10 use of just three such rotating gradiometers in fact provides a redundant 
measurement. This redundancy permits a least-squares best fit solution, 
further improving the quality of data obtained. Any DC offsets such as low 
frequency drift in flux gates or the fixed offset of a SQUID, for example, can 
also be determined and monitored to provide information about the operating 

15 conditions of the detectors. Such offsets can also be separated from gradient 
data in the Fourier domain. 

In preferred embodiments of the third and fourth aspects of the invention, 
data obtained is converted to the Fourier domain. In the Fourier domain the 
gradient information is distinguishable from the field information because the 

20 signals due to the gradient fields are centred at twice the rotation frequency 
while those due to the fields are centred at the rotation frequency. 

Further, information about the gxy component of the gradient tensor may 
be distinguished from that due the diagonal gradient components, even at the 
same frequency, because one involves taking the Fourier transform of a cosine 

25 waveform, which is real and even, while the other involves the Fourier 
transform of a sine wave, which is odd and imaginary. The Fourier 
transformation may be implemented by means of a conventional fast Fourier 
transform algorithm or even simple discrete implementation via a coarse 
sampling of say 8 data points per cycle of the rotation. 

30 To first-order any mismatches in the sensitivities of the field sensors or 

their relative alignment gives rise to contributions only at the fundamental 
frequency (discussed in more detail in the following) and therefore such 
mismatches do not determine the ultimate CMRR of the gradiometer. The final 
CMRR performance is determined by the resolution and precision of the data 

35 handling and the ultimate noise floor of the detectors. Accordingly, preferred 
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embodiments of the present invention may provide for a CMRR improved by 
several orders of magnitude over that of state of the art instruments. 

Further, since the external gradient and field are modulated by the 
rotation, but DC offsets, drift or very low frequency noise in the detectors are 

5 not, the rotation provided by the present invention facilitates separation of the 
these effects in the Fourier domain so that the DC and low frequency 
performance of the instrument is not limited by the performance of the 
detectors in these respects. 

Furthermore, in the case of magnetic gradiometry, if the gradiometer is to 

io be employed in a moving platform such as an aircraft, for example, a 
suspension system is generally used to isolate the instrument from vibrations 
which would otherwise induce significant noise by moving the instrument in the 
earth's magnetic field. However, in previous systems, such suspension 
systems have proven to be extremely difficult to engineer as significant 

15 isolation is required at sub 1 Hz frequencies. However, the rotating 
gradiometer provided by present invention modulates low frequency signals 
with the rotation frequency, thus facilitating the use of suspension systems 
which then only need to provide effective vibration isolation near the 
fundamental and first overtone of the rotation frequency. Such suspension 

20 systems are typically significantly easier to engineer than sub-1Hz suspension 
systems. 

In preferred embodiments of the third and fourth aspects of the invention, 
the at least three gradiometers are rotated at differing frequencies. Such 
embodiments of the invention facilitate separation of the frequency of operation 
25 of each gradiometer, and further facilitate a reduction in "cross-talk" between 
the respective gradiometers, of the type set out in the preceding discussion. As 
such cross-talk effects will be occurring at differing frequencies, the effect of 
one gradiometer operation on another can be filtered out or removed in the 
frequency domain. 

30 Accordingly, preferred embodiments of the present invention may yield 

significant advantages over previous field detection devices and methods in a 
number of ways. First, rotation of the gradiometer breaks the nexus between 
engineering tolerances (in positioning two sensors in precisely parallel planes) 
and the ultimate CMRR performance. Further, providing three rotating 

35 gradiometers in accordance with embodiments of the third and fourth aspects 
of the present invention yields measurements of all five components of the 
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gradient tensor and all three the components of the field, with redundancy in 
the data, all by using only 3 field detectors. Additionally, the rotation of the 
gradiometer breaks the dependence of the DC performance of the instrument 
on the DC performance of the detectors. Embodiments of the present invention 

5 also facilitate absolute value measurements of both the gradient and field 
components using SQUlDs. Embodiments of the present invention may also 
provide real time information of the condition of the SQUID in operation, for 
example as to whether or not the SQUID has remained in lock. As noted in the 
preceding discussions, previous SQUID devices have only been capable of 

io making relative field measurements, not absolute field measurements. 
Furthermore, embodiments of the present invention may improve the sensitivity 
of field detection by shifting the working point of the detectors to a higher 
frequency, where the noise performance of the detectors and the vibration 
isolation afforded by suspension systems is superior. Yet further, certain 

15 embodiments of the present invention can serve to reduce the effects of cross- 
talk by use of distinguishable rotation frequencies. 

Brief Description of the Drawings 

By way of example only, preferred embodiments of the invention will be 
20 described with reference to the accompanying drawings, in which: 

Figure 1a illustrates a low temperature superconducting gradiometer in 
accordance with the present invention; 

Figure 1b illustrates a high temperature superconducting gradiometer in 
accordance with the present invention; and 
25 Figure 2 illustrates a three-gradiometer magnetic field detection device in 

accordance with the present invention. 

Description of the Invention 

Figure 1a illustrates a low temperature superconducting gradiometer 10 
30 in accordance with the present invention. Gradiometer 10 comprises a first 
field vector sensor in the form of superconducting wire loop 11 Gradiometer 10 
further comprises a second field vector sensor in the form of superconducting 
wire loop 12. Sensors 11 and 12 are positioned substantially coaxially so as to 
sense field vectors along axis 13. . Further, sensors 11 and 12 are connected 
35 in a differencing arrangement, whereby a signal produced at terminals 14 
stems from the difference between the magnetic field vectors sensed by 
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sensors 11 and 12. In accordance with the present invention, gradiometer 10 
is configured to rotate about axis 15, which is at a normal to the page as 
viewed in Figure 1a. Rotation of the gradiometer is indicated by direction 
arrows 16. 

5 Figure 1b illustrates a high temperature superconducting gradiometer 20 

in accordance with the present invention. Gradiometer 20 comprises a first 
field vector sensor in the form of a SQUID 21. Gradiometer 10 further 
comprises a second field vector sensor in the form of a SQUID 22. SQUIDs 21 
and 22 are mounted on a suitable magnetically transparent mount. Sensors 21 

10 and 22 are positioned substantially coaxially so as to sense field vectors along 
axis 23. Further, sensors 21 and 22 are connected in a differencing 
arrangement (not shown), whereby a signal thus produced stems from the 
difference between the magnetic field vectors sensed by sensors 21 and 22. In 
accordance with the present invention, gradiometer 20 is configured to rotate 

15 about axis 25. 

Figure 2 illustrates a three-gradiometer magnetic field detection device 
30 in accordance with the present invention. The gradiometers are 
represented by cylinders 31, 32, 33. Cylinders 31, 32, 33 are respectively 
positioned on axes 41, 42, 43, and each of axes 41, 42 and 43 is not parallel 

20 with either of the other two axes. Each of gradiometers 31, 32, 33 is rotatable 
about unique non-parallel axes in accordance with the present invention. 
Operation of such a device is set out in more detail in the following. 



25 



30 
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.This calculation models the sensitivity of the gradiometer instrument to the magnetic field, and; 
gradient tensor in the presence of 6 pairs of misalignment angles introduced -by; tilting the;. / 
SQUIDs through arbitrary angles. The effects of precession are not considered explicitly but can 
probably be modelled to sufficient accuracy within the theory. It is demonstrated that the 
misalignment errors lead to contamination of the gradient and B-field components by terms that 
depend on the misalignments, but the effects of misalignment can be removed by the theory by 
introducing appropriate calibrations constants. A final solution is derived in terms of instrument 
coordinates and the transformations to geographic coordinates is also described. 

The resulting instrument provides a measurement of the true-value of the gradient , and field • 
components which is relatively robust with respect to changes in the calibration constants. 
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I. OVERVIEW 

The .present embodiment comprises of a set of three rotating axial gradiometers or drum..gra-\ 
diometers. In the 2-SQUID implementation the output voltages from a pair of SQUIDs operated'- 
in flux-locked loops are differenced to provide information about the magnetic field density and . 
the first-order magnetic gradient tensor. If the SQUIDs are perfectly aligned with identical sen- 
sitivities then the difference voltage provides no information about the magnetic field density but 
this information is restored if the sensitivities of the SQUIDs are mismatched by a known amount- 
The difference voltage also provides a measurement of the first-order axial gradient provided that \ 
the contributions from higher-order gradients are sufficiently small to be ignored. This is easily \ 
seen by expanding the B-field as a Taylor series about the origin which is take to be centre point 
of the axial gradiometer. 



A, Taylor Series for functions of three variables. 

The Taylor series for a scalar function B (x, y> z) about the origin is 

• b(w) _£^ +s | + ,|)"b ( o>. a) 

n=0 . x 

In the case of a 3-dimensional vector field the previous equation applies to each of the vector 
components, so the a-component of B at the point r, B a (r) , and at the origin, B a (0) , are related 
by 

B a (r) = B a (0) + (. §- x + y | + *!) B a H- ± + y£ + *|) B " 
If the second- and higher-order gradients of B are. small then 2 may be approximated by ■ 

- b °<°> + (*! + »! + *!K "" "* ; 

provided that r is small. 
Let 

8 



9ap 
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denote components of the the first-order gradients of B. Then eq. ?? becomes simply 

B a (r)=? B a (0) + xgax + y9cty + zg az (4) 
or in coordinate independent form 

B(r) = B(0)+?r . . (5) 

where g is the gradient matrix and r is the displacement vector to the point at which the B field 
is measured. 

Using 5 if ri and r are two points the difference in the B-field at these points is 
• B(n)-B(r ) = B(0) + ?ri-(B(0)+fr )=?(r' x -r') 

which is proportional to the product of the gradient matrix and the displacement between the 
points. 

II. INSTRUMENT AND DRUM COORDINATES. 

The drum gradiometers rotate about axes with orientations denoted by the vectors U, V and 
W in an umbrella configuration to be described in detail below. The axial gradiometers within the • 
drums are located at the end points of U, V and W, and facilitate absolute value measurements 
of linear combinations of some of the components of the magnetic gradient tensor and magnetic 
field with respect to a drum coordinate system, C', which is fixed to the axis of rotation of- each 
of the three drums. The systems of equations in drum coordinates are transformed to a common 
instrument coordinate system 3 C, and solved to yield absolute value measurements of the magnetic 
gradient tensor and magnetic field in instrument coordinates. 

For use in geophysical exploration the final measurements have to be expressed in terms of . 
the geographic coordinate system, Cq- This is easily done if the orientation of C with respect to 
Cq is known but in general this will require additional information from position and orientation, 
sensors. It is possible! /that, in some applications, the real-time measurements of the magnetic 
field provided by the instrument itself, might be used to provide information about the 

orientation of the instrument with respect to the magnetic field of the earth. Transformations of the 
final measurements between the instrument coordinates and geographic coordinates are described 
below in section VI. 
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A. Instrument coordinate system, C. 

The instrument coordinates system, C, is a right-handed Cartesian coordinate system defined 
as follows. The x-y plane is defined by the base of the instrument. The ^-direction is normal to 
the x-y plane and forms an axis of symmetry for the instrument in the sense that the three drum 
gradiometers are distributed symmetrically about the 2-axis/ Unit vectors in the directions of the . 
re-, y- and *-axes will be denoted i, j and k, respectively. 

As mentioned above the axial gradiometers are located at the end points of U, V and W and- 
rotate about the axes defined by these vectors. The direction of U is the vector obtained by.' 
revolving a unit vector i with its base at the origin of C about the —y axis through- an angle A.-. 
The directions of V and W are obtained by revolving U about the z-axis through angles of 120° •'" 
and 240° respectively. U, V and W then form a type of umbrella configuration being disposed 
symmetrically about the z-axis and tilted above the x-y plane through an umbrella afi^e* angle A." 
Hence 



U = (L u cos A, 0, L u sin A) (6) 
V = J -il^cos A, ^-L v cos X y L v sm A J (7) 



W 



= (^-~^L W cos A, ~^-L w cos A, L w sin A^ 



(8) 



in which L u , L v and L w denote the distances of the axial gradiometers from the origin of C. 
1. Properties o/U, V and W. 

• Note that these unit vectors form an orthogonal system only if 

U«V = V- W = w.u = o 

i.e. if 

7 - ^ cos 2A = 0 
4 4 

A = i arccos | = 0.615 48 rad = 35. 3° (9) 
2 ■ o 
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• In the instrument A = 45° so from equations 6, 7 and 8 

U=(cosJ,0,sin^)=.L u (V2 0 V5) ( 10 ) 
V- (-\cosl ^.cos£, sin^ -X.(-JV5 JV5 M) ■ (") 

W- (—omJ, -^cosj.sinj) =L tt (4V2 -JVS M) (12) 
2. Unit drum vectors 

It will prove useful to define the unit drum vectors, uo, vo and w 0l that are normal to U, V . 
and W respectively. Vector u 0 is the unit vector obtained by revolving a unit vector i with its. 
base at the origin of C about the +y axis through an angle (90 - A) degrees. Vectors v 0 and wo ' 
are obtained by revolving u 0 about the z-axis through angles of 120° and 240° respectively; In 
instrument coordinates 

uo= (sinA 0 -cosA) (* 3 ) . 

v 0 =(-IsinA ^(sinA) -cosa) ■(")'■ 
wo = ( -I sin A (sin A) -cos a) < 15 > 

N:B. 

u 0 -U=( s inA 0 -cosa) • L u (cos A, 0, sin A) =0 

v 0 -V=^_I s inA ^(sinA) -cos A ) • L v cos A, ^ cos A, sin A j =0 
wo • W= ( _I S i n A (sin A) - cos A ) • L w cos A, cos A, sin A^ =0 
as required. 

B. Drum coordinate systems C 

As mentioned earlier we introduce 3 drum coordinate systems, C", C" and C" associated with: 
the each of the three drums rotating about U, V and W respectively. Let x', y' and z' denote 
unit vectors in the directions of the x\y f and z' axes of the drum coordinate systems. C is defined 
so that x' = uo, y' = j, z' = V/L u and the origins of C\j and C coincide. With these definitions 
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C'u is a right-handed coordinate system which would become equivalent to C if the umbrella angle 
was set to 90°. 

C" and C" are then defined by rotating C about the *-axis through angles 120° and 240° so 
that in C", x" = v 0 and z" = V/L v while in C", x'" = w 0 and z'" = W/L w . 



1. Unit gradiometer vectors. 

The instantaneous direction of the rotating axial gradiometers will be defined, by unit gradiome- 
ter vectors u, v and w, which remain normal to U, V and W and undergo continuous rotation 
about the end points of U, V and W with angular frequencies u> u , oj v and w w respectively. 0 U , 
0 V and 6 W wiU denote respectively the angular displacements of the vectors u, v and w at time t,: 
from the positions of u 0 , v 0 and w 0 defined above, so that 0 U = u/ u t, 6 V = u) v t and 9 W = u w t. In 
terms of the drum coordinates C' v and G' w the unit gradiometer vectors u', v' and w'assuine. 
the same form 



u' = (cos 6 Ui sin 0 Ui 0) 
v' = (cos 6 Vl sin 0 Vl 0) 
w' = (cos 6 W , sin 9 Wi 0) 



(16) 
.(17). 
(18) 



C. Coordinate transformations between C and the C systems. 

In general if P and P' denote the components of a vector in the frames C and C respectively 
and if C is mapped to C" by a positive (right-handed) rotation through an angle k about an axis 
k (k = x,y,z) then 

P 7 = ^(k)P ( 19 )- 

where Rk (*) denotes the standard rotation matrices 

^ 1 0 0 ^ 

Rx (*) = 0 cos k sin k • ( 20 ) 

^0 -sin/c cos/cy . , ■ 

^cosk 0 -sin«\ 
Ry(K) = 0 10. . . ' . (4) 

^ sin k 0 cos k J 
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Jl, to- 



cos k sin k 0 \ 
— sin/c cosk 0 
0 0 \ } 



(22). 



Also if p and are the components of a Cartesian tensor in the frames C and G* respectively 
then 



9* = Rk{K)gR T k («). 
For example a rotation about the z-axis through angle k is 

9l 9* 9*z 



\9l3 0*3 033 / 



cosrc sinK 0 
— sin/c cosfc 0 
0 0 1 



( 9xi gi gia^ 1 

91 9 9 3 
\ 913 9 3 033 J 



cos/c — sin« 0 
sin k cos K 0 
0 0 1 



and hence 



9n = 1 (011 + 9 )-M(0il-0 ) cos + sin 2k 
0' = i (011 + 0 )- 1 (011-0 ) (cos 2k) -gi sin2« 
033 = 033 

01 = — (011 - 0 ) ( sin 2k ) + 01 cos 2k 
0 / 3 : = 0 3C0S«-pi3sinAc 
013 = 013COS/C + 0 3sin/c 



(33)' 



(24). 



1. Transformations C C. 

When — uo,y' = j and z' = U/L u , C is mapped onto C by a rotation about y through an 
angle (7r/2 - A) . Hence if B and B' denote the components of a vector in the frames C and C f 
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respectively then 



B' = 



B = 



sin A 0 —cos A 



0 1 0 
cos A 0 sin A J 

( sin A 0 cosA^ 

0 10 
y — cos A 0 sin A J 



B 



B' 



(25) 



(26)' 



N.B. If A = tt/4 then 



/IV2 0 -V*\ 
B' = 0 10 B (27);. 

while if A = 7r/2 25 and 26 give that P' = P, i.e. the coordinate systems C and C' v coincide, as 
required. Also, using equation 6, 



sin A 0 —cos A 



0 1 .0 
^ cos A 0 sin A J 

( sinA'O -cosA\ 



U 



0 10 

cos A 0 sin A 

. 0 

V 1 / 



L u cos A 



^ L u sin A j 



as required. 

Using 26 the unit gradiometer vector 



u 



/ COS0 U \ 

sin 6 U 

\ 0 J 



becomes 





^ sin A 0 cos A ^ 


f cos O 




u = 


0 10 


sin# u 






^ — cos A 0 sin A j 


V 0 ) 





sin 6 U 
— cos^tiC° s ^ y 



(28) 
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in C, Observe 



( cos B u sin A \ / L u cos A \ 



u-U = sin6 u • 0 =0 
y — cos 6 U cos A J 

' so u 1 U as required. 

a. Gradient Transformation, The transformations of the gradient tensor from G to C follower ■ 
from the use of the transformation matrix in 25, with the general transformation equation, 23, i.e. 



;• (29): 





A 


9'xA 




^sinA 0 


- cos A ^ 


^ 911 91 


913 ^ 


/ 


sinA 0 cosA^ 


A 


9* 


9*3 




0 1 


0 


91 9 


9 3 




0. 1 0 : 




9'* 


933 J 




^cosA 0 


sin A j 


\ 013 9 3 


933 J 


V 


-cos A 0 sin A j 



011 = 1 (011 + 033 - 011 cos 2A + 533 cos 2A - 2513 sin 2A) 
9' =9 

033 = - (011 +033 + 011 cos 2A - 533 cos 2 A + 2pi3 sin 2A) 



pi =V5foi -53) 
5' 13 = ~~ ( 2 013 cos 2A - 511 sin 2A + 533 sin 2A) 



In the special case A = 7r/4 



5 / 3 = 1 V2(5i +5 3) 
5u = -911 ~ 513 + -933 



9' =5 



533 = i 511 + 513 + A 533 
g[ = V2(01 -5 3) 
5i3 = ^511 - "533 
5 / 3 = 1 ^(5l +5 3) , 



> . 



(30). 



(31). 
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2. Transformations C *=* C". 



When x" = vo and z" = V/L v , C is mapped onto C" by a rotation about ^ through an angle . 
27r/3 followed by a rotation about y through an angle (it/ 2 — A) . 

^-^sinA ^v^CsinA) -cosA^ 

! -M - 1 o 

^--cosA -\/3(cosA) sin A j 

( — IsinA -i-v/3 -icosA \ 
!>/5(smA) -i iV3(cosA) 
- cos A 0 sin A 



B" = 



B = 



If A = w/4 



V 



B" = 



B" 



/ 



-JVS - 1 o 



B. 



Using 33 the unit gradiometer vector 



/ COS0„N 

sin5 v 
\ 0 J 



becomes 

( -^sinA -iyfe -±cosA \ 

1^5 (sin A) -± ^(cosA) 

y — cos A 0 sin A 

in C. 

Observe that 



(32)<: 



(33) 





^ cos 0 V ^ 


/ 




sin# v 




/ 


\ 0 / 


V 



f -^smXcos6 v - ^\/5(sin0 v ) ^ 
— - sin 0 V 4- -V5 (sin A cos 0 V ) 
— cos A cos 9 V 



(34). 



vV=0 



WO 2004/015435 



20 



PCT/AU2003/000999 



confirming that v 1 V as required. Also using 6 

/ — ^sinA ^VsCsinA) -cosA^ 

^ — icosA ^V^cosA) sin A J 



V 



= L V 



0 

V 1 / 



\ 


^ — ilr„COS A 


\ 




^L v cosA 




/ 


^ L v sinA 


) 



as required. 

a. Gradient Transformation. The transformations of the gradient tensor from C t&C follows, 
from the use of the transformation matrix in 32, with the general transformation equation, 23, i.e. 

\ ( -isinA -.*cpsA \ 

^(sinA) -1 ^(cosA) 



( g » r g » g» 3 \ /-isinA ^(sinA) -cosA^ / 



Jl Jit J" 

9\ 9 9 3 
V$13 9% 9h) 



-1 0 
-^cosA ^ (cos A) sin A 



511 91 913 
91 9 9 3 
\9l3 9 3 933 J 



— cos A 



" sin A . 
(35), 



J 
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011 = iflii + fa + X m - 3511 cos 2A - |fl cos 2A + ipsa cos 2A 
+i 513 sin2A-^i V5+ Jffl \/3cos2A-^ 3 V3sin2A 

g" - + isi a/3 + V3 (is! + Jffnv^) 

033 = iflll + |s + i S33 + 5 5ll cos 2A + | S cos 2A - ^33 cos 2A 
-ipi 3 sin2A-^i V3- 3 5l V3cos2A+ ±g 3 \/3sm2A 

(jj = 1 (cos A) g 3 - 1 (sin A) g\ + 1 (cos A) pi 3 \/3 
+|(sinA)pnV3-l(sinA)p >/5 

s « 3 = i 913 cos 2A + \gii sin 2A + |p sin 2A - *gn sin 2A 
— -g 3 \/3cos2A - \g\ \/3sin2A 



5 " 3 = -1 (cos A) pi - 1 (sin A) g 3 + \ (cos A) jmv/3 
-| (cos A) g V3- 1 (sinA)fli3\/3 



In the special case A = tt/4 



9" =1511 + 35 + A 5l V3 

533 = 1511 - A 513 + §5 + 1 533 - M51 + A \/35 3 

g>{ =-^51 +\^/2g 3 + lVQgu + \V&9i3-l^9 

913 = 8511 + §5 - I 533 - iVSffi 
5" 3 = -iV25l ~3\^5 3+iv^511-iV^5l3-iV65 , 



> . 



(3.7)*'. 
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When x'" = w 0 and z'" = W/L w , C is mapped onto C" by a rotation about z through an 
angle 4?r/3 followed by a rotation about y through an angle (ir/2 - A) . Hence 

''-AsinA -i-\/3(sinA) -cosA^ 

1^/3 -1 0 B 

^ —- cos A -->/3 (cos A) sin A J 
( -IsinA ±y/3 -*cosA ^ 

-1^3 (sinA) -i -iV3(cosA) B'" 
^ -cosA 0 sinA J 



B'" = 



B = 



(38).;. 



(39). 



If A = ir/4 



B"' = 



-i Q 



B. 



(40) 



Using 39 the unit gradiometer vector 



w"' = 



becomes 



/ _I s inA iy/S -±cos\ 



w = 



\ 



V 



-±V3(sinA) - 1 -iv^(cosA) 
-cosA 0 sinA J 



cos 0,0 
sin^ 

V 0 / 

''cosO 
sin0 v 

V 0 / 



/ -± sin A cos 0 V + i^sin -*v) ^ 
— 1 sin 6 V - 1 V3 (sin A cqs 0 v ) 
— cos A cos &v J 



(41) 



in C. 

Observe that 



wW=0 
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confirming that wlWas required. Also using equation 8 





( 


- sin A 


1 fn /_!_ \ \ 

— ~V3 (sin A) 


— cos A 






w'" = 




M 






n 
u 


W 






l- 


- cos A 


\6 \COS A) 


sin A j 








/_ 


-sin A 


-i\/3(sin A) 


— cos A ^ 




/ 


—-L w cos A ^ 










0 






-^L^cos A 




\- 


1 cos A 


-A (cos A) 


sin A j 




\ 


L w sin A J 






(°) 














0 





























as required. 

a. Gradient . Transformation. The transformations of the gradient tensor from ' G to • G UL fol- . 
lows from the use of the transformation matrix in 38, with the general transformation equation;, 
23, i.e. 

( --sin A 

& -1 0 

y -I cos A -^(cosA) sin A y 



011 9i 013 x 

./// _/// jut 

01 0 9 3 

\ 013 0 3 033 / 



^ ^ 011 Pi 013 \ / - 2 sinA -4 cos A \ 



91 9 9 3 
\ 513 5 3 533 / 



-^(sinA) -1 -^(cosA) 
^ -cos A 0 sin A J 



(42) 
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011 = 5511 + |s + i S33-i0ilcos2A-|$ cos 2A +^33 cos 2 A 
+^l 3 sin2A+ \/3-^i \/3cos2A+^ 3\/3sin2A 

^=B^i + l^ + i 533+ipiicos2A+|p cos 2A -^33 cos 2 A 
-ipi 3 sin2A+^l V3 + ^i V^cos2A- i p 3 \/5siii2A 

g'{' = i(cosA)p 3 -i(sinA)pi - 1 (cos A) pm/3 
-^(sinA)pnV3+^(sinA)p n/3 

= ip 13 cos 2 A + sin 2A + |p sin 2A 
-lp 33 sin2A+ A s 3 V3cos2A + ^i v / 3sin2A 

p" 3 = -i (cos A) 5i - 1 (sin A) 5 3 - ^ (cos A) ffn\/3 
+\ (cos A) 5 V3+ ±(sinA)<>i3\/3 

In the special case A = tt/4 

011 = iffll + i »l3 + |ff + 1 <733 + |\/3ffi + i V3ff 3 

^33 = 5511-^13 + 15 '+^33 + ^51 - A V3ff3 
g»' = -|V2pi + 4%/2 S 3 - iv%l - |%/6ffl3 + iV6ff 
ffl3=iffll + 55 - ^33 + ^1 

p" 3 = - 3 - i^Pll + i%/6ffl3 + , 



> . 



(43) ; 



(44) 



D. SQUID Configuration 



The SQUIDs composing the axial gradiometer rotating about U axe located at the endpoints. 
of the vectors r = U ± pu, where p = d/2 and d is the total distance between the SQUIDs. In 
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drum coordinates 









^pcos 




r' + = U' + pu' = 


0 


+. 


p sin 0 U 










I 0 J 


\ 




f °1 




^ p cos 0 U ^ 


( 


r'_ = U'-/ro' = 


0 




p sin 0 U 






\L U J 




I 0 j 


\ 



p COS 9 U ^ 

psxn6 u 



Using eqs. 45 and 46 in equation 5 the magnetic field at the SQUIDs may be written 

(BL\ f g' xx jL g' xz \ (pco S e u \ 



B'+ 04) - 



'a; 



psin0 u 



%c 9yy 9yz 
V fix Sy* 9zz ) 
f B' x + ludxz + P ( sin *«) fl^ + P (cos 0 U ) Si« ^ 

B' y + L u g'yz + P ( C0S 0u) 9'xy + P ( sin Q u) 9yy 

\B' z + p (cos 0„) g' xz + p (sin 0 U ) ff^ + / 



B'_ (r'_) = 



\B' Z J 



( sL &y 9' x z\ (-PCOSBU^ 



— psin0« 
L u J 



9yx 9yy 9yz 
\ 9'zx 9zy 9'zz J \ 

( B' x + L u g' xz - p (sin 9 U ) g l xy -p (cos 6 U ) g'^ > 

B'y + ~ P (COS 6 U ) g' X y ~ ? (SUl 6 U ) g'yy 

^B' z -p (cos 0„) ^ - p (sin 9 U ) g' yz + I-u0 zz ) 



in which 



\B' Z J 



is the magnetic field at the origin of C 



(45) 



(46).. 



(47)' 



(48). 



(49) 



E. SQUID outputs. 

In the instrument the gradiometer signal is produced by forming the sum 

U = U++U- 



(50) 
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where C/+ and E/_ are the responses of the SQUIDs when operatedin a flux-locked loop. In response 
to the magnetic fields in equations 47 and 48 (7+ and U~ may be written 

U- = S f - B'_ (r'_) + a- . (52) ; 

where the Si and S' denote the SQUID sensitivity vectors. These are normal to the SQUIDs with : 
lengths proportional to the sensitivity of the SQUID system (V/T) and the a are arbitrary offset 
voltages which can change value upon each reset of the SQUID system electronics. In a perfectly 
aligned system Si = Siu' and S' = -S u' but in any real system the SQUID Vectors will be tilted-' . 
with respect to u' due to alignment errors. To characterize the alignment we introduce 2 pairs of 
alignment error angles. Let cr\ (<r ) denote the error angle between the projection of the Si (S 7 )•• 
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vector onto the x'-y' plane and the unit vectors u' (-u') , and r\ (r ) is the angle of inclination of 
Si (S' ) above the x'-y f plane. Then 

f Si COS T\ COS 9\ \ . 

cos rising , ^i = (^ + ^i) (53). 
Si sin n J 



SU" 



\ 



/ -S COST costf ^ 
-S cos r sin 6 , 0 — (#u + o ) 
£ sinr ^ 

If cri = <r = n = r =0 then = S\v! and S 7 = -5 u' as required. 
■ To simplify these introduce 

Mi ^ Si cosri, iVi = Si sin n 
M ~S cost , N = S sinr 



(54) 



(55) 



so that 53 and 54 become 



S' = 



t Mi cos $i * 
Mi sin 6i 

, * J 

( ~M cosO ^ 
— M sin# 



(56) 



(57). 



Using 47 and 56 in equation 51 gives 

/ Micosfli^ ( 
Misin^x 
Ni ) 



U+ = 



\ 



V 



B' x + L u g'xz + P ( sin e u) g' X y + P (cos 0 U ) gL 
B' y + L u g' yz + p (cos 9 U ) g' xy + p (sin 0 U ) gyy 
yB' z + p (cos 6 U ) g' xz + p (sin 8 U ) cf yz + Lug' zz J 

= (Mi cos Ox) (B' x + Lj xz + p (sin 9 U ) g' xy + p (cos 6 U ) g' xx ) 

4 (Mi sin Ox) {B' y + L u g' yz + p (cos 9 U ) g'^ + p (sin 9 U ) g' w ) 

+ (JVi) (B' z + p (cos 0 U ) g' xz + p (sin 6 U ) g' yz + L u g' zz ) + a+ 



+ a+ 



(58). 

(59) 
(60) 
(61). ; . 



= B^Mi cos 6x + L u g' xz Mx cos 6x + g'xxPMx (cos 0„) (cos 0i) + g'^pMx (sin 6 U ) (cos 0i) . 

+ B' y Mx sin «i + ^£uMi sin <?i + g' w pMx (sin 0 U ) (sin 0i) + ^Mi(cos0 u ) (sin 9x\ . 

+ B' z Nx + g zz NxL u + g xz NxP (cos e u )+g yz Nxp(sm6 u )+a + ■ . (62)u 
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where we have used that fey = Qyz, 9xz = 9zx and g yz = Pzy from our knowledge of the properties 
of the magnetic gradient tensor. 

Similarly using 48 and 57 in equation 52 gives 



/ 


— M cos 9 


) 


( 




-M sine 






\ 


N 


) 





B f x + L u g xz - p (sin 0 U ) g'xy-p (cos 0 U ) 0** 
B[, + Lugy Z - P (cos 6 U ) g f xy - p (sin 6 U ) g 1 ^ 



+ 0- 



= (-M cos 9 ) (B' x + L U 0^ - P ( sin e «) 9 xy - P ( cos *u) 9 f xx) 
+ (-M sin 0 ) (B^ + Lupyz - P ( cos $4y - P (sin 0„) 
+ (N ) (B' z - p (cos 0 U ) g' xz - p (sin 9 U ) g' yz + L u g zz ) + a_ 



(63) " 

(64) , 
.(65) 

(66). 



= -B' X M cos 0 - L u g' xz M cos Q 4- M pp^ (cos 0„) (cos 0 ) + M pg xy (sin 0 tt ) (cos 6) 
— M By (sin 0 ) - g' yz M L u (sin 9 ) 4- $yi/ M p (sin 0 tt ) (sin 0 ) 4- g'^M p (cos 0 tt ) (sin 0 ) 
4- B^JV + <4iV L u - JV p (cos 0 U ) - fl^JV p (sin 0 tt ) + a_ ' (67) 
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Substituting 62 and 67 into 50 gives the gradiometer output 

U = B X M\ cos 0i + Lug f xz Mi cos 0i + g xx pMi ( cos ft*) (cos 0i) + g^pMi (sin Q u ) (cos 0\\ 

- B X M cos 0 - L u g f xz M cos 0 + M pgt^ (sin 0 U ) (cos 0 ) + M pg xx (cos 0 U ) (cos 0 .) 
+ B^Mi sin $i + Qy Z L u Mi sin 0i + S^pMi (sin 0 U ) (sin 0i) 4- 0^Mi ( cos ft*) ( sin 

— M By (sin 0 ) - g' yz M L u (sin 0 ) 4- c/^M p (cos 0*) (sin 0 ) + g'yyM p (sin 0 tt ) (sin 0 ) 
+ B Z N! + g' zz NiL u + t/^Wip (cos 0 U ) 4- c^ip (sin 0 U ) 4- a + 

+ 4- c/^iV I u - p (cos 0 U ) - g' yz N p (sin 0 U ) + a_ > : . : >' : . 

= B' x Mi cos (0 U 4- <ri) 4- Lug' xz Mi cos (0 U 4- <r x ) + P^Mi (cos 0 W ) (cos (0 U .+ cri)) . ' 
+ g'sypMi (sin 0 U ) (cos (6 U 4- cri)) - cos (0 W 4- a ) - L u g' xz M cos (0 U + cr ) 

4- M pc^, (sin 0 U ) (cos (0 U + a )) + M pg'^ (cos 0 U ) (cos (0 U + * )) 
+ JB^Mi sin (0 U + ai) + s^L u Mi sin + a l) + s'yyPMi (sin 0 U ) (sin (0 U 4- <ri)') 
+ p^pMi (cos e u ) (sin (0 U + tn)) -M B f y (sin (0 U + a )) 

- p^M L u (sin {Q u 4- cr )) + t/^Af p (cos 0 U ) (sin (0 tt + cr )) 

4- g'yyM p (sin 0 U ) (sin (0 U 4- cr )) 4- B^Ni 4* g l zz NiL u + S^lP (cos 0«) 
+ <7v*^iP (sin 0 U ) 4- a+.+ B^JV 4- s^JV L u 

- P'( cos ft*) - Py*^ /° ( sin ft*) + a- 

= B^Ma (cos cri cos 0 U - sin cri sin 0 U ) + L u ^Mi (cos a\ cos 0 U - sin a\ sin 0 U ) 

4- ^Lc^^l ( cos °"1 cos ft* - sin cri cos 0 U sin 0 U ) + g^pM\ (cos cri cos 0 U sin 0 U - sin cri sin 0 U ) 

— B' X M (cos cr cos 0 U - sin cr sin 0 U ) - L u g xz M (cos a cos 0 U - sin cr sin 0^) 

4- M pg f xy (cos cr cos 0 U sin 0 U -sin cr sin 0 U ) + M pc;^ (cos cr cos 0 U - sin cr cos0 u sin0 u ) 

4- B f y Mi (cos cri sin 0 U 4- sin cri cos 0 U ) + g yz L u Mi (cos cri sin 0 U 4- sin cri cos 0 U ) 

4- g f yyPMi (sin cri cos 0 U sin 0 U 4- cos ai sin 0 U ) 4- piypMi (cos cri cos 0 U sin 0 U 4- sin cri cos 0u) - 

— M B y (cos cr sin 0 U 4- sin cr cos Q u ) - g' yz M L u (cos cr sin 0 U 4- sin cr cos 0 U ) 

4-p^M p(coscr cos0 u sin0 u 4-sincr cos 0 U ) + p^M p (sin cr cos 9 U sin 0 U 4- cos cr sin 0 U ) 
4- B^iNTi + p^^iL u 4- g' xz Nip (cos 0 U ) + g f yz Nip (sin 0„) 4- a + 
4- B' Z N + p^iV L u - <4iV p (cos 0 U ) - g yz N p (sin 0 U ) + 
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= B X M\ cos ai cos 9 U - B' x Mi sin a\ sin 9 U + L u g xz M\ cos ai cos 9 U 

- L u g xz M\ sin c\ sin 0 U 4- g' xx pM\ cos <ji cos 0 U - g xx pMx sin cri cos 0 U sin 6 U 
+ g'xypMi cos ci cos 0 U sin 0 U - g^pMi sin cri sin 6 U -B X M cos a cos6 u 

+ B X M sine sin 9 u — L u g xz M cos a cos0 u 4- Lyjg xz M sincr sin0 u 

+ P9xy c ^s c cos0 u sin#u - M pp^sino- sin 9 U + M pg xx cos ar cos 0 U 

- M P9xx sin a cos sin ®u 4- -ByMi cos c\ sin 0 U + B^Mi sin a\ cos 0 U 

+ g' yz L u M\ cos crx sin 0 U 4- g' yz L u M\ sin ai cos 6 U 4- g'yypMi sin <7i cos 0 V sin #u 
+ ^/>Micos<7isin ^ u + p^pAficosaicos^ u sin^ u + ^pAfisincricbs 9 U 
-M B f y cos o- sin 9 U — M B f y sin cr cos '9 U - Ir U cos a sin 0 U 
-~ £ u sin<7 cos0 u 4-p^M pcos cr cos 9 U sin 0 W + 'g'^M p sin a cos 0 U 

/?sincr cos0 u sin0 w 4-0^,M pcoscr sin 0 U 
+ B' z Ni 4- pi^i A* 4- g' xz N iP (cos 0 W ) 4- ^Wi/? (sin 9 U ) 4- a + 
4- B£ W 4- L u - g f xz N p (cos - g' yz N p (sin 0 U ) + a_ 
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= B X M\ cos 6 U cos cri - B X M cos 9 U cos cr 
+ B X M sin cr sin 0 U - B X M\ sin cri sin Q u 
+ jB^Mi cos o"i sin 6 U - B^M cos a sin 0 U 
4- By Af i sin cri cos 0^ - B f y M sin cr cos 6 U 
4- I'u^-^i cos cri cos 0 U - LuP^-^ cos a cos 0 U 
^uP^z^l sin 0*1 sin 0 U + Lug xz M sin cr sin# u 
+ Pyz-^u^l cos cri sin 6 U - g f yz M L u cos (7 sin 0 U 
+ g^^u-Mi sin <7i cos 0 U — <7yz £u sin a cos 6 U 
+ B^iVi + B^N + p^iiu + g' zz N L u + a+ + a, 
+ P^xP^l cos (7i cos 0 U + M pg'zn cos C7 cos 0^ 
-J- g'yypMi cos cri sin 9 U + p^M p cos (7 sin 0 U 
+ g f xy pM\ sin (7i cos #u + gicyM p sin a cos 9 U 

— g xy pM\ sin (7i sin 0 U — M pg^ sin (7 sin 0 U 

— p^pMi sin cri cos 0 U sin Q u — M' pgt^ sin cr cos 0 U sin 0 U 
4- g'yypMi sin cri cos 0 U sin 0 U 4- g'yyM p sin cr cos 6 U sin 0 U 

4- 2g' xy pMi cos cri cos #u sin 6 U 4- 2M cos cr cos 0 U sin 6 U 
4- g f xz N\p (cos 0 W ) - c^iV /o (cos 0 U ) 
+ 9yzNlP (sin flu) - g' yz N P (sin ^u) 
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= (Micoscri — M cos cr )B x cobO u 
4- (M sin cr - Mi sin 01) B x sin9 u 
4- (Mi cos <r\ — M coscr )B f y s\xi9 u 
4- (Mi sin cri M sincr )B y cos9 u 
+ X u (Mi cos (ri — M cos cr ) cos 0 U 

— Z w (Misincri — M sin cr ) sin 0 U 
4- jL u (Mi cos cri — M coscr )p^sin^ u 
4- L u (Mi sin ai - M sin a ) p^, cos 9 U 
+ L u (Ni + N)g' zz 

+ B' Z (N! + N) 
4-a 

+ £ (Mi cos 0*1 4- M cos cr ) g' x (cos 0 U ) 
4- p (Mi cos cri 4- M cos cr ) (sin 0 U ) 
4- p (Mi sin ai 4- M sin cr ) g xy (cos 0 U ) 

— /0 (Mi sin cri 4- M sin a ) g^ (sin 0 U ) 

— p (Mi sin cri + M sin a ) <^ (sin 9 U ) (cos 0 U ) 
4- p (Mi sin cri -h M sin cr ) ^ (sin 9 U ) (cos 0 U ) 
4- 2p (Mi cos cri 4* M cos cr ) g 1 ^ (sin 0 U ) (cos 0 U ) 
+ P9' XZ (cos 9 U ) (iVi-iV) 

+ ^(sin^)(iVi--^) 

where a = a+ 4- a_. Expanding the trigonometric functions within each terra and collecting 
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harmonics gives 

U= (L u g' zz + B' z ) (N x + N) + a 
+ ip(Micoscri + M cos a )(g' xx + g' yy ) 

•f (Mi cos cri - M cos a ) (i?^ + LttP^) cos 8 U 
+ (Mi sin (7i - M sin <7 ) (B^, + Lutf^) cos 
+ p(N 1 -N)g' xz cos0 u 

- (Mi sin <7i - M sin cr ) {B x + L u g xz ) sin 0 U 
+ (Mi cos <ri - M cos cr ) (jB^, + sin 0 U 

+ p(Ni -N)g , yz 8m0 u 

+ |p (Mi cos (7i + M cos cr ) (p^ - cos 20 u 
•f p (Mi sin <7i + M sin a ) p£. y cos 20 u 

+ \p (Mi sin ax + M sin cr ) (g^ - p^) sin 20 u 

+ p (Mi cos cri + M cos cr ) g'^ sin 20 u . (68) 
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Hence 



i7 = Mu (LuQzz + B^) + a 
+ f u (pL + ^iy) 

+ (-Si + L u g' xz ) cos 0 U 
+ *?U {By + Lu9y Z ) cos 
+ «ui4rCOS0u 

+ <n(B , y + L u g' yz ) sin0 u 

+ £u (pja - %/) cos 20 u 
+ 2(Xtt)^cos2«„ 

+ X«(Pw-5i»)sin2^u 

+ 2(f 1t )fl^8in2tf tt . (69) 



where 



£ u 4 ^p(Mi cos (7i + Af cos cr ) 
Xu — -p(Misincri -f M sin cr ) 



Cu — (Ml COS (71 — M COS (7 ) 

rj u = (Afi sincri — M sin <7 ) 



(70). 



N.B. In the absence of alignment errors, i.e. a\- a = 0, n = t =0 from equation 55 



Mi = Si, JVi = 0 
M =S , N = 0 
and from equations 70 Xu> *?u, v u and ^vanish, and 




[71). 
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Finally if Si = S = S say, then <; u also vanishes and the only non-zero constant' is 

& = " (*?)■ 

III. DATA EXTRACTION 

One method to extract the data is to take a series of measurements over a- number of angular 
increments of the pick-up loop structure and solve the resulting system of equations algebraically. 
A series of voltage measurements is now taken at angles 8 U = — 7r/4, 0, 7r/4, 7r/2, 37r/4> 7r, 57r/4,. 
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37r/2 to yield the following set of equations: 

Uq Au(0=: -tt/4) = a + B' 3 fMu + 0' n f« + g' n Xu - 2g[ U + 0' £u - ^ Xu + ^33/^ 

C/ AU(6 = tt/4) = a + B^ u + ^„ - g' n Xu + 2^ & + ^ & + + £u033A*u " 

U A ±U{d = 3tt/4) = a + B^ u + p' u C u + s'nXu - 2pi £ u + s' £u - Xu + i^Mu 
- A iWl3 W2 " 1 W3W2 - i ^u5l3^V / 2 +' 1 

I/ fl , AU(0 = 5?r/4) = a + B^/iu 4- « - SuXu + 2p'i & + 0 ; f u + 0' Xu + ittfSs/** . 

+iBi7 ?u V5-^7 ?u% /2-^u%/2- 1 B'c uV ^- Vi3*W5- V 3 <W5 ( 7 3) 

Wi3%V2 - 1 W 3 - 1 Wl3^V2 - iWaCuVS 
C7 X 4 27 (0 « 0) = a + B' ^ + Bfa + B£/i u + 2g f n ^ u + ffi 3 v„ + 2pi X u 

U Z ±U(6 = tt/2) = a - «Bfr u + B'<r u + B^u - 2^ + 0' 3 «u + 2p' 
-^vPl3^u + Lug'tfu + £u033Mu 

C/ 5 4 C/ (0 = tt) = a - B' ?7u - Bfc tt + B^u + 20n£u - 9i3 v u + Vi Xu 
—Lug'tfu - L u ^' 13 Cu + ^uP33Mu 

U 7 AU(6= 3?r/2) = a + B'^u - B'<r u + BsA*« " Xu - S>u + 2^ £u 
+-&tipi3*7u-^uff'3^i + -^tt^3Mu 
Note that even with all of these misalignment errors considered the following L.C. vanishes if 
the SQUID stays in lock throughout the measurement series: 

U0-U1 + U -U 3 + U4~Us + U 6 -U 7 = 0 / ( 7 - 4 >.. 

So this is still a good check to test the SQUID lock and the integrity of the data with respect to.! 
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random pointing errors which might occur during the measurement sequence. 



A. Solving for gradient components. 

Assuming that £ u , Xu> Cu, Vu> ^u, f^u are all known by calibration, the 8 equations in 73 contain 
10 unknown independent variables, namely B x1 B f yy B Z) g xxi g^ g zzi g'^ g yzi g^ a, since .voltages 
Uo,„Uj are determined by measurement. The set of equations in 73 is clearly- uriderdetefmihed 
and additional data from the other two drums is necessary to determine a complete ; solution. " 
Nevertheless it is possible to solve for the off-diagonal gradient component, g xy , and the difference . 
between the diagonal components, p^. — g'yy, independently of the measurements from the other 
two drums. This is done by forming the following linear combinations of eqs. 73 to eliminate the 
dependence on the SQUID offset voltage a and on the field components, 



(75) 
(76) 



where 



Then 



provided 



fa' N 



A I A J J 
~9xx~ 9yy- 

/ Xu(Ul-U3 + Us-U7)-tu(Uo-U + U4-Ue) \ 
8(£„ + Xu) 

Xu(Uo - U + Ui - Pe) + gu (Pi -U3 + U5- Ui) 



(e«+x«)^o 



(77) 



(78) 



which is always satisfied because both £ u and Xu are real and £ u > 0. 

The results from these additional gradiometers give solutions with exactly the same form in 
their respective coordinate systems, i.e. 



8(£„ + Xv) 
4«„ + Xv) 



(79)' 
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and 



v A "7 



/ Xw (Wi -W3 + W5- W 7 ) - £ w (Wq - W + Wi -W 9 ) '\ 
&ttv> + Xw) 



(80) 



Xv>(W 0 -W +W i -W 6 )+t w (W 1 -W 3 + Ws-W 7 ) 
*«• + %.) 

Equations 77, 79 and 80 are now transformed to the instrument coordinate system and- the. 
resulting system of equations may be solved to yield the gradient measurements. 



1. Transformations for the "proof of concept" instrument. 

In the case of the present embodiment A = 7r/4, and the transformations for the xy^x and 
yy gradient components of the C", C" and C" coordinate systems are those given by' the; 'relevant' 
eauations in 36, 37 and 44, which are brought together here for convenience: 

ffi = -|\/2ffi + \ V2g 3 + § v^ffn + 3^13 - IVeg 
iff = -\y/2gi + IV^g 3 - gv/6511 - \V^913 + iV&9 



(81). 



ffil ~ ~9\\ — 913 + -933 
9' =9 

9ii = gsu + l 9t3 + iff + 1 933 ~lV3gi -'V3>3 
9" = fffn + \g + 1 gi V5 

511 = Sffll + 1 913 + Iff +^33 + 1^1 + 1 V / 5ff3 

ff"'=lffll+^ff -iv^ffl 

Using the last six of these equations we can derive the transformation equations for the diagonal 
component differences A', A" and A"', i.e. 

A' = ffW 



= Tjff" ~ 913 ~ 9 + 2^33 

X a 1 A 
= -A a --A -<m 



(82) 



A if _ Ji n 
= 911-9 

5 1 1 3 yr 1 /r 1 

= - J Ai 4- ^ A3 - 701 V5 - J<7 sV3 + ^13 



(83). 
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= 011-0 
5 

"8"" ' 8 



5 1.1 3 /- 1 n I 
-Q0H +o9+ o033 + V3 + -g 3 V3 + -513 



= - jAi + |a 3 + ^01^ + ip 3 >/3 + i^ia 



(84) 



where 



(85) 



Ai =gn-9 
A =0 - 533 

A3 - 533 - 511 

Combining these with the first three equation of 81 yields the following set . of equations 

9'l =V2 5 l -V253 
g'{ = -\s/2gi + \y/2g 3 + | VUga + |>/6Ai 
5?' = -|V2pi + |V2 fl 3 - \\fcgi3 - |\/6Ai 
A'=iAi-iA -513 

A" = -|Ai + ^A 3 - {y/Zgi -^VSg» + 1 ga 
A'" = -|A X + U 3 + |V3ffi + ^5 3 + *5l3 . 

Although the corresponding matrix of coefficients, with respect to the variables 
91 ,9 3,913, Ai, A , A 3) 







0 


0 0 


0 > 




(91 \ 




f-5i \ 








0 


0 




9 3 




-5i 








|%/6 0 


0 




513 




"5i 


0 


0 


1 


_l 1 


0 




Ai 




"A'x 






_1 


i ° 


_i 




A 




-A' 


\-lv3 


-1V3 


_1 


i 0 










I -A3 J 



is singular a solution can by obtained after observing that, by definition, 



Ai + A + A3 = gn - g +0 - 033 + 033 - 011 = 0 
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and use this equation to augment 86, i.e, 

9 >» = -\^/2gx + \s/2g 3 - - |V6Ai 

A'-±Ai-±A -<m 

A" = ~\ Ai + ±A 3 - f V3<?i - 1 V3 5 3 + l 9\z 
A'" = -£Ai + ±A 3 + |\/3pi 4- A Vfy 3 + 1 513 
0 = Ai + A + A 3 

The corresponding matrix form 

/-±V2 ±V5 0 0 0 n \ 

JV5 -JvS -\y/i o 

0 0 1 

cr 1 

8 



3^ 1^/3 _1 1 

i 

8 



-1 



1V3 



o o 



i 

0 
0 

-1 -I) 



fa, \ 



91 
9 3 
913 
Ai 
A 

VA 3 y 



~9\ 

J" 
~9i 

-A' 

-A" 
-A'" 

0 / 



is transformed to a square matrix by multiplying on the left by the transpose of the coefficient 
matrix to give 



(f 


3 


0 


0 


0 


o\ 




(91 \ 




f Vi - \g'{ V2 - | A" V3 - Jtff ^ + |A'"V3 > 


3 


9 
1 


0 


0 


0 


0 




9 3 




- Vi + K - 1 A" V3 + Jflf V? + -A"'n/3 


0 


0 


I 


-J 


I 


1 




513 




-a + i a" + aa"' + is? v/e - ya V6 . .;• 


0 


0 


-I 


f 7 


3 
4 


I 




Ai 




-A' - 1 A" - \ A'" + }rf V6 - jj»f V6 


0 


0 


1 


! 


5 
4 


1 




A 




-iA' 


v° 


0 


1 


I 


1 






\A 3 J 




^ I a" + lA"' . / 



which is now non-singular. The inverse of the coefficient matrix is 



/ 


\ 


-« 


0 


0 


0 


o A 




16 


44 
75 


0 


0 


0 


0 




0 


0 


136 






6 

""5 




0 


0 


88 


sot 




44 




0 


0 






436 
~ 5 


68 

-75 


\ 


0 


0 






-* 


hi) 
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(87) 



and the solution of the system of equations, 86, is 

/ (A'" - A") + -if (2g\ - g'{ - sf) ^ 

-$(A"/-A") + £ (g'{ -25i + 5 f) 
I5 (A"-2A' + A'") + f (g'{ -gf) 

TgA' - & A" - £ A'" + ^ vs - a* V6 

-ft A - ^A" - &A'" - „rf V6 + ntf V6 
V ? A' + |A" + fA'"- 15 ffi'V6 + I5 Pi"V6 / 

Note that at this point this solution has been obtained without using the tracelessness of -. the 
gradient tensor and is written in terms of the differences between the. diagonal components. -The.;:; 
diagonal components, however, follow immediately from equations 85 and the traceless .of the 
gradient tensor, i.e. given that 



9 3 
913 
Ai 
A 
\A 3 / 



911-9 = Ai 
9 - 933 = A 
£33 ~ 9ll = A 3 
911 + 9 +933 = 0 , 



1 -1 0 

0 1 -1 
-1 0 1 

1 1 1 



9 

\933 J 



A 
A 3 

\ 0 



(88) 



then 



( gn ) 




^-A 3 + Ai^ 


9 


_ 1 
~ 3 


-Ai+A 


\933 J 




^ -A + A 3 J 



(89). 



Substituting for the component differences from the last three equations in 87 gives 

/ \ / _^ (_A' + 8A" + 8A'") 4- fof - 9*{' ) ^ 

(11 A' + 2A» + 2A'") 4- (<?f - <?? ) ( 9 °) 
^ §(A' + A"-f A'") • / 

The final solution for the gradient components follow after substituting for the A's and./s; 
from equations 77, 79 and 80 into equation 90 and the first three equations of $7. The resulting 



911 
9 

\933 J 
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equations are brought together here for convenience. 

911 = -45 ("A' + 8A" + 8A'") + ^ (g'i - gf) 

g = (11A' + 2A" + 2A'") + -4 (<#' - g'{ ) 

S33 = i(A' + A" + A"') 

91 (A'" " A") + (Vi " 9'{ ~ 9'i! ) 

j 3 = ^ (A'" - A") + £ (g'i - 2g[ + gf ) 

giz = is (A" - 2A' + A'") + 4 (g'i - g'i') 



in which 

, _ Xu(Ui -U3 + U5- U 7 ) - tu(U 0 - U + Ui - U 6 ) 
91 " 8(€ U +Xu) 



„ _Xv(Vi-V3 + V 5 -V 7 )-Z v (Vo-V +V 4 -V e ) (91) 
51 _ 8(^ + Xv) 

,„ _ Xw(Wi-W 3 + W s -W 7 )-U(Wo-W + Wj-W 6 ) 
91 ~ 8(^ + xJ 

., Xu(Uo-U +Ui-U6) + UUi-U3 + Us- u t) 
A_ 4(^ + Xu ) 

.„ xv(Vo-v +V4-V6) + V3 + V5-W) 

_ X^O-W +W4-W6) + g«(Wi-W3 + W B -W r 7) 
A ~ 4(^ + xJ 

and 

£ u 4 ip (Mi cos (71 + M cos cr ) 
Xu 4 ±p (Mi sin cri + M sin C7 ) 
Mi 4 Si cos 7i, Ni = Si sin n 
M =S cost , N = S sinr 
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£, Special case: "Proof of concept instrument 11 with zero alignment errors. 
If all of the alignment errors vanish then 

Xu = 0 

and equation 91 reduces to 

SU = "35 (-A' + 8A" + 8A'") + -£ (g'{ - g'{> ) 



g =- 2g (iiA' + 2A" + 2A"') + -^( ? i"-«7i') 
S33 = |(A' + A" + A'") 

gi = -@(A'»-A«) + -£ (2fli -s? -<#') 
53=^ (A'" - A") + £ ^ - 2 ff i + 5 f) 
013 = 15 (A" - 2A' + A'") + ^ (g'{ - ) 
and 

, -(Uq-U +U 4 - U 6 ) 
91 = 8fu 

-(V Q -V +V 4 -V 6 ) 
91 = 



,„ - (Wo - W +Wi-W 6 ) 
91 = 8& 



4£„ 



4^ 



. . . (92). 
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IV, SOLVING FOR FIELD COMPONENTS. 

If the gradiometer is intentionally unbalanced, i.e. 5i ^ S , then it is capable of returning' 
information about the magnetic field in addition to the gradient. As for the gradient case considered 
above we form linear combinations of 73 to both eliminate the dependence on the SQUID offset 
voltage, a, and to reduce the number of gradient terms. To eliminate the SQUID offset we start 
with the differences between voltages corresponding to opposing positions, and then take the sums: 
and differences of these pairs, i.e. 

(Uo - U 4 ) + (U - U 6 ) = 2rj u B' V2 + 2<; U B / 1 V2 + 2v u g[ 3 \/2 + 2L u 77 u p' 3 y/2 + 2£ u Cusi 3 \/2 

(Uo - C/4) - (U - Uq) = 2r} u B' 1 V2 - 2<; u B l y/2 - 2v u g l ^ + 2L u Wi3V^ - 2L u <r u / 3 v / 2 ■ \: " 

(Ui - U 5 ) -f (U 3 - Ur) = 2B'r} u - 2B' 1 Tj u 4- 2Bi<r u + 2B f <; u 

+ 2^13^ + 2p' 3 f u — 2L UJ 95[37?tx *+• 2L u g' 3 Tj u + 2£ u £i 3 c u 4- 2£ u £' 3 c u 

(^1 - - (^3 - ^7) = 2B , l r ]u + 2J3'ifo 4- 2JBi«u - 25' c u 

+ 2pi 3 u u - 2g' 3 v u 4- 2L u g[^7] u 4- 2L u g r 3 rj u 4- 21^^^ — 2L u ^' 3 c u 

and then simplify further by forming the following L.C. of these equations: 
Ui + Us-Us-U 7 Uq + U -U4-Ue 



2 2V2 
Ui-U z -Us+U 7 U 0 -U -U 4 + U 6 



= Cu^' - 77u^i + ^'3 ~ LuVu9l3 + LySujg's 
= t?u5' 4- <r u £i 4- v u gi3 + W3 + ^uCuPi 3 ' (94) 



2 2V2 
As for the gradient calculation previously the results from the other two drums have exactly the 
same form in their respective coordinate systems, i.e. 

Vl ± V3 _ V5 _ V7 _ Vo ± V -v 4 - v 6 m ^ _ + vvg% _ Wg + ' { ^ 

yi _v 3 V5 + v 7 _ v 0 ~v-v 4 + v 6 = + + ^ + w , + £wftA • 

2 2V2 

tVl+^3-W 5 -W 7 TVp + W -W A -W 6 „, ,„ 

2 2V2 = *° ~ VW 1 

+ v w g r 3 — L w T] w gi3 + £u>Su#" 3 (97) 
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2 2^5 1 

Equations 93 - ?? may be simplified by making the following substitutions: 

{/T = — ^3 (L u <; u + V U ) + Wl 3 7?u 



-(98). 



2 2\/2 



2 2\/2 

2 ~~ 2y/2 

y X a Vi ~ ^3 - ^5 + V? _ Vo-V -T4 + V6 

w% A Wi- W3-W5 + W 7 W 0 -W -W± + W 6 



+ Lv9\$nv — g%v v — L v g n tf; v 
- 9isv v - L v g%7j v - IrvP^Cv 
+ L w gi^rj w — g%v w — Lwg%Sw 
— 3 u iu — Lw9 f/ 3Vw — £u>0i3Su> 



(99) 



2 2v^ 
Observe that all of the terms on the right hand sides of these equations are known i.e. the U, V 
and W voltages are returned by instrument, the gradient components are known in the instrument 
frame and can be easily found using coordinare transformations, and the remaining terms are 
known by calibration. With these substitutions 93 - ?? become 

= <; U B' -rjuB^ 

W* = <; w B'"-7i w B f {' 
W* = r ho B"' + <; w B'{' j 

Equations 100, corresponding to measurements in the C", C" and C M coordinate systems, may 
be solved, by transforming to the common instrument coordinate system using the coordinate-, 
transformations for field components derived earlier, but the resulting sets of equations rapidly 



(100) 
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become intractible. It is simpler to first solve these pairwise in tlieir respective coordinate systems,. 

1 fu^ u -U^u\ 



(101) 



i / vk v -v*w 



(102); 




(103) 



and then carry out the transformation to the instrument system. As derived above the relevant 
transformation equations for the field components are 

^ sin A 0 —cos A ^ 
B'= 0 1 0 B 

^ cos A 0 sin A j 
B[ = - (cos A) B 3 + (sin A) Bi 

B' = B • • - 

B£ = (cosA)Bi + (sinA)B 3 

^-^sinA -\/3(sinA) -cosA^ 
B"= -1 0B 

^— ^cos.A ~\/3(cosA) sin A y 
' Bi / =-(cosA)5 3 - 1 (sinA)jBi + i(sinA)B V3 

< B" = -*B - i^iv^ 

^ B'I = (cos A) Bi + (sin A) B 3 + 1 (cos A)B \/3 
f -I gin A -i>/3(sinA) -cosA^ 
B'" 1 V3 - 1 0B 

^-icosA -iVS(cosA) sin A y 
' Bi" = - (cos A) B 3 - 1 (sin A) B x - 1 (sin A) B \/3 

< B'" = --B + i Biv / 3 
Bf = -± (cos A) Bi 4- (sin A) B 3 - ± (cos A) B \/3 
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and if A = 7r/4 



B' =B 

Bf- JV5(-Bi + >/5S -2B 3 ) 
B" = -±(BiV3 + B ) 
■Bf = (-Si - %/3B -2B 3 ) 
B'" = i (V3Bi - B ) 

and equations 101 - 103 become 



= iy/2(B 1 -B 3 ) 



= B 



= \y/2(-Bi + y/lB -2B 3 ) 



_ 1 



(BlV3 + S ) 



= JV5 (-Bl - V5B - 2B 3 ) 



(1.04) 



(105) 



A least-squares best-fit solution may be found by first rewriting this overdetermined system of 
equations in matrix form in terms of the unknowns B\,B , B3 to give 



( l sft 


0 






Jlu + ?u. 
U*T) U + c/ku 


0 

-JV5 


1 

_1 


0 ! 
0 






.»7u + S 

>+s 








\B 3 ) 








_1 


0 y 




Jlw + ?tu. 
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and then multiplying on the left by the transpose of the coefficient matrix to yield 

2 Vv + *v . 2 »7w + ^ 



4 Vw + ^w • 2 + 2 Vw + Sw 
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We bring the final forms of the solutions here for convenience 

Uku ~ U^u Vg - V tffr 2y/3 + V^v) 



B _ 4 y*tfa + yt^ | ^6 y*^ - 2 y + y 

9 Vu + Su 9. ?7„ + ^ 9 Tfo-fCt, 
9 ^ + <; w 9 Tfo,+ ^ 



where 



(108) 



yT = —j= 4- L v g l3 Tj v - g 3 v v - L v g 3 Cu 

V+ = —7= gi$v v - L v g S 7] v - L v g l3 <; v 

^ 4 Wl + Wa - W5 - W7 _ w^w-Wt-We + ^ _ ^ _ ww 

^ 4 _ W 0 -W 2 -W< + W 6 _ g%vw _ Lwgll , vw _ ^ 

. sis = x m - 1 933, <?" 3 = - i^<? 3 + IVegu - {Vegia -\^g . 

{/a = l V2(9i +93), S13 = &n + fs - i S33 + i\/3</i 

si's = |su + fa - A S33 - $V5gi , <?" 3 = -| V201 - \^g 3 - i >/5*u + |V6<m +. jy4 . 

Formally the scaling coefficients, v Ul <; u and % are defined by equations 55 and .70 but would 

normally be determined by calibration. 
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V. EXAMPLES 

A. Precession about the z-axis. 

Say cri = a = 0, Si = S = 5. If the gradiometer precesses about the the axis of rotation by 
virtue of an error angle of 0 say, then the theory developed above suffices to model this by letting^ 
t*i = — t =<f>. 1 . Using equations ?? 



£u — pQS cos (j) 1 

v u — 2pQSsin(j) 
Ck = 0 
7?u = 0 



a u 4 (2Q 5 /? cos <£) 1 

A* = (O)" 1 
tft" 4 0b - + E/4 - ^6 
f/* 4 t/i - Us + C/s - U 7 
Mi^Scos^, TVi -^sin^ 
M ^Scosfa N = -£sm<£ 



(110) 



VI. APPENDIX 1: GEOGRAPHIC COORDINATE SYSTEM, C G . 

The geographical coordinates system, Cqi is defined as follows: 
#<3-axis: The x^-axis is along the traverse in the direction of increasing station number ("grid north"). 

yc?-axis: The 2/G-axis is at 90° to the rccr-axis, parallel to the ground and points right as you look' 
along the traverse ("grid east"). 

ZQ-Bxis: The zoaxis is vertically down (pointing at the centre of mass of the earth). 

This is a right-handed Cartesian coordinate system in terms of which a unit vector with azimuth 
A (the angle measured clockwise from grid north) and inclination I (angle measured positive 
downward from the horizontal) is 

unit vector with azimuth A and inclination I = (cos A cos J, sin A cos /, sin J) . (HI) 



1 If the gradiometer drum is tilted through and an angle <j> in the a>z plane say, the pick-up loop vectors are also 
tilted through the same angle and the positions of the pick-up loops are altered. The change in position has a 
negligable effect on measurements of gradients produced by distance sources. Since by definition t is the angle of 
the sensitivity vector above the x-y plane, a tilt though <j> implies that n = <f> and r = -<f>. 



WO 2004/015435 



51 



PCT/AU2003/000999 



1. Transformations between C and Cg- 

The transformation between C and Gq assumes a particularly simple form if the instrument is 
oriented with its base horizontal and its x-axis aligned due north. In this case C is mapped to Cq 
by a positive rotation through an angle 7r radians about the z-axis and the transformations of the 
field and tensor and data become simply 



9G = R^gRn 1 



where 



Hence 



R* = R* 1 — 



and 



B G 
\BG3 J 



f 9G11 9G1 9G13^ 

9G 1 9G 9G 3 
\ 9G31 9GS 9G33 ) 



1 0 0 
0-10 
0 0-1 

/ B a \ 
—B 

/ 911 -91 

-9 19 
\ -931 93 



-913^ 
9 3 
933 ) 



(112) 
(113) 



.(114)' 



(115) 



(116) 



More generally if the instrument has azimuth A and inclination J (in geographic coordinates) 
then C is mapped to Cq by first rotating about y through an angle of -I then by a rotation about 
the resulting z through an angle A followed lastly by a rotation about the resulting x of 7r: Hence 



B G = Rx (ir) Rz {A) Ry (- 1) B 



( Bqi\ 

b g 

\BG3 J 



\ 



( 1 0 0 

0 COS 7T sin7T 
0 -sin7T COS7T 

f cos A cos I sin A cos A sin J 

sin A cos I — cos A sin A sin J 

y sinJ 0 -cos/ J 



\ ( cos A sin A oW 
-sin A cos A 0 
0 0 1; 

b 

\B3) 



cos(-J) 0 -sin(-J) 
0 1 0 



(117) 



\ 


(bA 




B . 


/ 


\Bij 



(118) 
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and 



/ 0(711 9G1 0G13^ 

9G1 9G 9G 3 
\ 0G13 9G 3 0G33 / 

( cos A cos I sin A cos A sin 7 \ / 
sin A cos 7 —cos A sin A sin J 
sin J 0 -cos 7 



011 91 013 

01 9 9 3 

\ 913 9 3 933 J 
In the special case A = 7 = 0 we get so that 



\ 



\ ( cosAcos7 sinAcos7 sin7 \ 
sin A - cos A . 0 
cos A sin 7 sin A sin 7 -cos 7 



Bi 
-B 
\-B 3 / 



V 



\ 



9G1X 9G1 9G13 
9G1 9G 9G 3 
\ 9G13 9G 3 9G33 J 



-91 -013 ^ 



011 

-01 9 9 3 
\ -013 0 3 033 J 



as above. 
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The approximation used in obtaining equation (3) above, which is 
significantly simplified as compared to equation (2), involves truncating the 
Taylor series expansion of the magnetic field as a function of position to include 
the effects of the B-field and the first-order gradients only. This simplifying 

5 approximation is desirable because the Taylor series expansion has an infinite 
number of terms. The approximation is justified because the magnitudes of the 
higher-order terms decreases as the order of the terms increases, for three 
reasons. The higher-order terms in the Taylor series are scaled by prefactors 
which may be written as the product of two factors. Firstly, one of these factors 

10 is a scaling term which decreases with increasing order. The second factor is 
proportional to 1/r", where r is half the distance of separation between the field , 
detectors, and n is the order of the term. Since the separation between the 
field sensors is small, raising this to an increasing power rapidly decreases the 
magnitude of this factor. The third reason the approximation is justified is that 

15 higher-order gradients fall off faster with distance from the field source than 
lower-order gradients, so the signal at the gradiometer is dominated by the field 
and the first-order gradient components. 

The effects of misalignments of the field detectors on the output voltages 
from the instrument are also scaled by these higher-order terms, but these 

20 effects are small and similarly decrease with increasing order. 

Referring to equations (6)-(8), it is noted that the value of A, may be any' 
value in the range of 0° to 90°, inclusive. As noted in equation (9), when X is 
approximately 35.3°, an orthogonal system is provided. In the presently 
described embodiment, X is 45°. 

25 It will be appreciated by persons skilled in the art that numerous 

variations and/or modifications may be made to the invention as shown in the 
specific embodiments without departing from the spirit or scope of the invention 
as broadly described. In particular, the present invention may be applied in 
conjunction with the invention set out in the co-pending PCT application 

30 claiming priority from Australian Provisional Patent Application No. 2002950624 
and filed simultaneously with the present application, the contents of which are 
herein incorporated by reference. The present embodiments are, therefore, to 
be considered in all respects as illustrative and not restrictive. 



